Mode sum regularization approach for the self-force in black hole spacetime 
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We present a method for calculating the self-force (the "radiation reaction force") acting on a 
charged particle moving in a strong field orbit in black hole spacetime. In this approach, one first 
calculates the contribution to the self-force due to each multipole mode of the particle's field. Then, 
the sum over modes is evaluated, subject to a certain regularization procedure. Here we develop 
this regularization procedure for a scalar charge on a Schwarzschild background, and present the 
results of its implementation for radial trajectories (not necessarily geodesic). 
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A promising source for the LISA space-based gravita- 
tional wave detector is the capture of a small compact ob- 
ject by a supermassive black hole. The waveforms in this 
scenario are expected to be well approximated by a model 
of a small pointlike mass (a "particle" ) moving near the 
black hole M . One also expects such a model to provide 
valuable information in more general scenarios involving 
strongly gravitating two-body systems. To accurately 
predict the orbital evolution of a particle near a black 
hole (thus allowing the design of templates needed for de- 
tection and interpretation of gravitational wave forms), 
it is necessary to include nongeodesic effects caused by 
the interaction of the particle with its own field. Namely, 
one must tackle the problem of calculating the self-force 
(sometimes dubbed the "radiation reaction" force) expe- 
rienced by a particle in a curved spacetime. The usual 
approach to this problem is via a perturbation analysis, 
in which the charge q of the particle (which may represent 
its mass or its electric charge, or — in the toy model con- 
sidered here — its scalar charge) is assumed to be much 
smaller than the mass M of the black hole, and one looks 
for the leading order (oc q 2 ) self-force effect. 

There has been extensive effort, particularly over the 
last few years, to develop methods for calculating the 
self- force in curved backgrounds. The main challenge in- 
volved in this calculation is to deal with the various infi- 
nite quantities associated with the divergence of the field 
at the particle's location, which necessitates the intro- 
duction of an appropriate regularization technique. One 
calculation method, pioneered long ago by DeWitt and 
Brehmc Q], deduces the particle's equation of motion 
by imposing local energy-momentum conservation on a 
world-tube surrounding the particle's worldlinc An- 
other approach j| is based on computing the flux to infin- 
ity (and across the horizon) of conserved quantities which 
are tightly related to the orbit's constants of motion, par- 
ticularly the particle's energy E and azimuthal angular 
momentum L z in the Schwarzschild and Kerr cases. This 
approach, however, is inadequate for analyzing the evo- 
lution of the third constant of motion in Kerr, the Carter 
constant Q. Also, such techniques do not account for the 
non-dissipative part of the self-force p,p). So far, cal- 



culations of the self-force in actual examples have been 
restricted to very few simple cases, such as static charges 
in the Schwarzschild and Reissner-Nordstrom geometries 
A general formal framework for obtaining equa- 
tions of motion for a test particle in a curved spacetime 
has been formulated recently by Quinn and Wald 
however, the practical implementation of this approach 
in actual calculations remained, so far, a challenging task. 

Previously, Ori proposed that radiation reaction 
effects may be calculated by evaluating the contribution 
of each Fourier-multipole mode £mu> of the retarded field 
to the radiative evolution, through the local force applied 
on the particle by this mode, and then summing over all 
modes. This approach has two advantages: First, each 
individual mode of the field turns out to be continuous 
(and the corresponding self-force to be bounded) even at 
the particle's location. Secondly, calculating each £mu> 
mode of the field becomes a relatively simple task, as 
it only requires the solution of an ordinary differential 
equation (DE). In |l(J this sum-over-modes approach has 
been proposed for the calculation of the adiabatic, orbit- 
integrated, evolution rate of the three constants of mo- 
tion in Kerr, i.e., E, L Zl and It might have been 
hoped that the same method would also be useful for 
calculating the momentary self-force; however, this naive 
procedure for calculating the self-force turns out inappli- 
cable: Although each mode yields a finite contribution, 
the sum over all modes is found to diverge. This is the 
situation even in the simple case of a static scalar charge 
in flat space: For such a charge, located at a distance ro 
from the origin of coordinates (with respect to which the 
spherical harmonics are defined) , the contribution of each 
multipole £ to the radial component of the self- force is the 
same: —q 2 /(2r^). Obviously, the sum over £ diverges. To 
overcome this type of divergence, one must introduce a 
certain regularization procedure into the calculation. 



"The mode sum for the adiabatic, orbit-integrated, evolution 
rate of E, L z was shown to converge pc| . It is not clear yet 
whether the corresponding mode sum for Q converges or not. 



1 



In this paper we describe the basic elements of a new 
regularization scheme for the mode sum. As this method 
does not involve any weak-field or slow-motion approxi- 
mations, it should allow effective calculations of the self- 
force for strong field orbits. Here, we outline the method 
for a scalar particle moving on a Schwarzschild back- 
ground, and present the final results for the special case 
of radial motion (not necessarily geodesic). We also men- 
tion some results, derived recently by one of us pi] , for 
the case of a circular orbit. A more detailed account, 
including a more systematic presentation of the regular- 
ization procedure, full details of the calculations involved 
(for radial trajectories), and generalization to a wider 
class of static spherically symmetric black hole space- 
times (including, for example, the Reissner-Nordstrom 
and Schwarzschild-de Sitter spacetimes), will be given in 
a forthcoming paper [fl2| . 

Schematically, the calculation of the self-force in our 
approach contains two (essentially independent) parts. 
In the first part, which deals with the contribution of the 
individual modes, the appropriate ordinary DE for each 
mode £mui of the field is solved (numerically, in most 
cases), and the contribution F^ u of each mode to the 
self-force is evaluated. (Alternatively, one may numeri- 
cally solve the 1+1 partial DE in the time domain, for 
each £ and m.) In the second part, a certain regulariza- 
tion procedure is applied to the mode sum. This part 
involves the analytic calculation of certain regularization 
parameters, using a local perturbative analysis. This 
manuscript deals with the second part, namely, the ana- 
lytic regularization scheme. The first part — the numeri- 
cal determination of F% mu — was implemented by Burko 
for several scenarios Jl3|— |l5|| - 

Throughout this paper we use relativistic units 
(with G = c = 1), metric signature ( — h++), and 
Schwarzschild coordinates t, r, 6, (p. 

The total self-force on a scalar particle is composed of 
three parts p,^6[: 



p(total) _ ci(ALD) i p(Rii 



p(ta.\\) ^ 



(1) 



The first term is an Abraham-Lorentz-Dirac (ALD)-like 
term, F a ALD ^ = ^q 2 (d a ~a 2 u a ), where u a and a a are, re- 
spectively, the four-velocity and four-acceleration of the 
particle, a 2 = apa@ , and an overdot represents covari- 
ant differentiation with respect to the particle's proper 
time r. The second term is a Ricci-related term, which 
vanishes in the Schwarzschild case considered here. The 
third term, F a tml \ which represents a non-local contri- 
bution to the self-force, may be expressed as 



F< tail > = lim Fi £ ) 



(2) 



where]]] 



t Following ]9j, we use here = g<ji iM as the basic expression 
for the force applied on a charged particle by the scalar field. 



G, a [x"]x$(rj\ dr. 
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Here, x^(t) denotes the particle's worldline, and the inte- 
grand is evaluated at x^ = x$ = x^ (t = 0) (the self- force 
evaluation point). G is the retarded Green's function, 
satisfying 



UGix^-x^) 



-47T 



5 4 {x» 



(4) 



subject to the causal initial condition, i.e., G = when- 
ever x^ lies outside the future light cone of x" 1 . [In Eq. 
([|), □ represents the covariant D'Alembertian operator.] 
From the practical point of view, the challenging task is 
the calculation of the tail part, which requires knowledge 
of the Green's function everywhere along the particle's 
past worldline. In what follows we therefore focus on the 
calculation of the tail term Fa ■ 

We decompose F a ** and G(x^]x'^) into spherical har- 
monic £-modes, 



F, 



(e) _ 



(5) 



(Fa and G l are the quantities resulting from summa- 
tion over the azimuthal number m). We can then write 

: r. - SF5 e) ) , (6) 



F<"» = Urn lim + E( J 



where Fa , SF a , and F^ denote the force associated 
with the £ multipole of the field sourced, respectively, by 
the interval r < — e, the interval r > — e, and the entire 
worldline (through the ^-mode Green's function G ). The 
quantity F^ may be identified with the sum over m, ui of 
the contributions from all modes £muj for a given £ (recall 
that in calculating the total field's mode £mui one takes 
the source term to be the entire worldline). 

It should be emphasized that the limit e — > and the 
sum over £ in Eq. (H) should not be interchanged (other- 

wise the entire contribution from SF a , which is crucial 
for the regularization procedure, would be lost). Also, a 
clarification is required here concerning the meaning of 
the last equality in Eq. (^): The quantity F^ is well 
defined at r = ro (the particle's location). The situa- 
tion with F^ and SFa is more delicate, however. Each 
of these quantities has well defined values at the limits 
r — > Tq and r — > Tq . However, these two one-sided values 
are not the same. Equation (|6|) should thus be viewed as 
an equation for either the limit r — > of the relevant 

quantities (i.e., F^ and SF a ), or the limit r — > of 
these quantities. Obviously, this equation is also valid for 
the averaged force, i.e., the average of these two one-sided 
values. Throughout this paper we shall always consider 
the averaged force. (Of course, the final result of the cal- 



culation, F a tal1 ^ , having a well defined value at the eval- 
uation point, should be the same regardless of whether 
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it is derived from its one-sided limit r — > Tq , or from 
r — > , or from their average.) 

Next, we seek a simple e-independent function h a , such 
that the series YlA^a ~ ^1) would converge. Once such 
a function is found, Eq. (o) becomes 

where 

D a ^ Km VfeW-d (8) 

In principle, /i^, can be found by investigating the asymp- 
totic behavior of as £ — > oo. It is also possible, how- 
ever, to derive h l a from the large-^ asymptotic behavior 

of SFa^ [the latter and must have the same large-^ 
singular asymptotic behavior (for fixed e), because their 
difference yields a convergent sum over £]. Obviously, 

in order to determine h e a and D a from 5Fa , one only 
needs the asymptotic behavior of SFa^ in the immediate 
neighborhood of e = 0. This allows one to derive h l a and 
D a using local analytic methods, as we now describe. 

To proceed, we define the (Eddington-Finkelstein) null 
coordinates, v = t + r* and u = t — r*, where r* = 
r + 2Mln(r/2M - 1). The Green's function at for a 
source point at x" 1 may be expressed as Jl2[ 

G l = LP e (cos X ) g ^ V ' U; f' U>) 6fa - - u'), (9) 

rr 

where L = £ + 1/2, O is the standard step function, Pi 
is the Legendre polynomial, and 

cosx = cos 9 cos 9' + sin0sin#'cos(<p — <p'). (10) 

The "reduced Green's function" g*(v, u; v' , v!) satisfies 

9%u + V e (r)g e = 0, (11) 

with the supplementary characteristic initial conditions 
g l (v = v 1 ) — g e (u = u') = 1. The effective potential 
is V e (r) = (f/4)[£(£ + l)/r 2 + f/r], where / = 1 - 
2M/r and a prime denotes d/dr. We re-express it as 
V e (r) = L 2 V {r) + Vi(r), where V = f/Ar 2 and V x = 
(//16r 2 )(4r/'-l). 

To explore the asymptotic behavior of the function g l 
for small spacetime intervals (and for large £), we next ap- 
ply the following perturbation analysis. We first Taylor- 
expand V (r) in the small deviation of r from tq. It is 
convenient to take the small expansion parameter to be 
r* — r*o- We also define 

A r = 2V^ 2 L{r r - r* ), (12) 

where Vbo = Vb(ro)- Expressing r* — r»o in terms of 
A r /L, substituting in the Taylor expansion, and collect- 
ing terms of the same powers in L (with fixed A r ), one 
finds 



V\r) = Vbo [L 2 + L (/iA r ) + (/ 2 + f 3 A 2 r )] 



+OQ./L), (13) 
where /i, f%, and /3 are coefficients given by 

/ 1 ^iv - 3/2 y = r 1 / 2 (r/'-2/), (14) 

h = V^Vi = r/' - 1/4, (15) 



h = ^V ~ 2 V = — [(ff /f + f") + 3(/ - rf). (16) 

Here, an overbar denotes d/dr*, and all quantities are 
evaluated at r = r . We now expand g l in the form 

oo 

g t = J2L- k g k (A r ,A rl> z), (17) 

fc=0 
1/2 

where A r / = 2Vq L(r* — r*o) and 

z = 2L[V 00 (v-v')(u-u')] 1/2 . (18) 

(The variable z represents the geodesic distance, scaled 
by L/tq, between the Green's function evaluation and 
source points, to leading order in r — r' .) Substituting 
the expansions (|l^) and ( |l7| ) in Eq. ( |i"l| ) and comparing 
powers of L, we obtain a hierarchy of equations for the 
various functions gk, of the form 

gk, V x + gk = S k - (19) 

Here, y = Vqq 2 L(v — v'), x = VqJ 2 L(u — u'), and Sk is a 
source term determined by gk><k- In the analysis below 
we shall only need the terms k = 0, 1, 2. For these values 
of k, the source terms are Sq = 0, Si = —fiA r go, and 
S 2 = -/iA r9l - (/ 2 + £^)So- 

The solution to Eq. ( |19|) for fc = 0, 1, 2, subject to the 
initial conditions gk{v = v') = gk{u = u') = 5ko (which 
conform with the original initial conditions for g l ), is 
given by [Q 

5o - Jo(z), (20) 
g 1 = --f 1 zJ l (z)(A r +A r ,), (21) 

.92 = -\zJi{z) [h{A 2 r + A r A r , + A 2 .,) + 3/ 2 ] 
+ ^z 2 J 2 (z) [3/ 2 (A r + A r , f - 8/ 3 ] 

+^/iVJ 3 (*), (22) 

where J„(z) are the Bessel functions of the first kind. 
Combining Eqs. (|^) and ([l7|), one obtains an expression 

for G e in powers of 1/L. One next constructs by 

,0+ 

5Fi^ = q 2 G e a [x»;x£(T)]dT (23) 
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(evaluated at = Xq). 

From this point on, we shall restrict attention to ra- 
dial trajectories. In this case, the only non-vanishing 
components of the self-force are F r and F t [note that 
Pi(cosx) = 1]- Expanding the various r-dependent 
quantities involved in the calculation (such as r' , A r /, and 
z) in powers of r, replacing r by L _1 (Lr), and collecting 
powers of X/L (for fixed Lr), one obtains an expression 

for 5Fa in the form of an expansion in X/L, 



5F i(e) 



LHW+Hj»+Hj?>/L + 0(L- 



dX, 



(24) 



(k) 

where A = — Lr. The coefficients H a are functions of A 
(but not of L or t otherwise), which also depend on the 
parameters ro, u r , if , and u r (evaluated at x^). These 
functions, all being certain linear combinations of a few 
terms A'J^A) (with integers i and n), will be given ex- 
plicitly in Ref. @. 

Considering the asymptotic form of Eq. (p3) for large 
£, one finds that the regularization function br a takes the 
form 

h f a = A a L + B a + CJL, (25) 

where[] 

A a = Hi°\X)dX 7 B a = H a 1 \X)dX, 
Jo Jo 

POO 

C a = / H a 2 \X)dX. (26) 
Jo 



Substituting Eqs. (|24j) and (|25j) in Eq. (g), one obtains 
D a = - lim V / LffW + iJ« + ffW/L dA (27) 



[the contribution from the 0(L~ 2 ) term in Eq. ( pi[ ) van- 
ishes at the limit e -> Q]. 

In conclusion, the tail part of the self-force is given by 

F (tail) = ^ (pt _ AqL _ Ba _ C J L j _ ^ (2g) 



motion, an expression of this form is also valid for any 
other trajectory. (Of course, the details of the trajec- 
tory under consideration will enter the calculation of the 
regularization parameters.) 

We shall now present the results obtained for the values 
of the above regularization parameters in the case of a 
radial motion (not necessarily geodesic). First, one finds 
|l2| that both parameters A a and C a vanish identically: 



A a — C a — 0. 



(29) 



In fact, it can be shown that at least the r component 
A r vanishes for all orbits in the Schwarzschild geome- 
try p| , and preliminary results indicate that the same 
holds for C r . (Recall that we are considering here the av- 
eraged force: the "one-sided" values of A a are, in general, 
nonzero [jl2j.) The vanishing of C a plays a crucial role 
in the regularization scheme: as it turns out ]l2| |, the pa- 
rameter D a cannot be well-defined [specifically, the limit 
e — > in Eq. (||) diverges] unless C a vanishes. Thus, 
the result ( p9|) may be viewed as a demonstration of the 
scheme's self-consistency. 

On the other hand, the parameters B a and D a 
are generically nonvanishing. For a radial motion in 
Schwarzschild, a calculation based on Eq. (|2^) yields JlJ] 

^(radial) = _ (( jr + ^ _ (3Q) 

zr 

The calculation of D a is more complicated. Basically, it 
involves transforming in Eq. ( |27j ) from a summation over 
£ to an integration over L f[ followed by an evaluation 
of the resulting double integral in the limit e — ► 0. This 
calculation, whose details will be presented in Ref. Eg) , 
yields, in the case of radial motion, 



-q 2 (a a - a 2 u D 



(31) 



Remarkably, D a is exactly the same as the ALD-like term 
F a AL . Thus, the contribution of D a to the tail term 
cancels out the local term in the expression for the total 
self- force, Eq. (|l]), leading to the simple result 



(32) 



(From the above construction of h a it follows that the 
sum over £ converges like ~ £~ 1 .) The implementation of 
our regularization scheme thus amounts to analytically 
determining the regularization parameters A a , B a , C a , 
and D a , using Eqs. ( p6| ) and (p7|). We point out that, 
although Eq. (pq) has been developed here for a radial 



Moreover, our analysis of a more general class of non- 
vacuum, static, spherically symmetric spacetimes (for 



which, in general, F a 1 ^ 0) yields G3| 
D a = F^ ALD) + i^( Ricci ) ; 



(33) 



* The integrals below include oscillatory terms which do not 
converge as A — > oo. In our calculation scheme we simply 
drop these terms. This will be justified in Ref. . 



§When performing this transformation, one must also take 
into account the 0(e 2 ) contribution arising from the difference 
between the sum over £ and the corresponding integral. J13] 
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implying that Eq. (|32j) is valid in this more general non- 
vacuum class as well. 

The regularization method described above has also 
been implemented for the case of a scalar particle in a 
uniform circular motion around a Schwarzschild black 
hole Here we merely present the results, concerning 
the r-component of the regularization parameters in this 
case. For circular orbits, the parameter D r (as well as A r 
and C r ) is found to vanish. |**| The parameter B r takes 
the form 11 



^(circular) 



2^7 



2I a - 



r - 3M 



ro 



2M 



(34) 



Here, 7 = (1 — V 2 )~ x / 2 , V is the tangential ve- 
locity with respect to a static observer, i.e., V — 
r a{— dtt)^ 1 ^ 2 (dtp I dt) (for an equatorial orbit), and I a = 
F(1/2,1/2;1;V 2 ) and I b = F(l/2,3/2; 1;V 2 ) are hyper- 
geometric functions. 

A special case of both Eqs. ( pp| ) and (pi) is the one 
corresponding to a static particle, for which we find 



r>(static) ^ 

r 2rg 



ro 



M 



ro 



2M 



(35) 



Recently, Burko numerically calculated F lmul for static 
p3[ and circular [jl4j orbits, and used the above regu- 
larization scheme to calculate the self-force. The above 
expressions for B r [Eqs. ( |34| ) and (p5j)], as well as the 
vanishing of A r and C r , are in excellent agreement with 
the large-^ limit of his numerically-deduced F l (obtained 
by summing F u over m and lj). 

To conclude, in this paper we have presented a practi- 
cal prescription for the calculation of the self- force acting 
on a scalar particle in the Schwarzschild spacetime. This 
prescription involves the (essentially straightforward) nu- 
merical calculation of the mode contributions F^, and 
the analytic derivation of four regularization parameters. 
The tail part of the self- force is then given by Eq. (28). 
For radial trajectories, the total self-force takes the sim- 
ple form ([32]) , with B a given by Eq. (^0|) . Based also 
on other cases studied so far one may, perhaps, 

conjecture that Eq. ( [32] ) holds for any trajectory in any 
spherically symmetric spacetime. 

The generalization of this mode-sum regularization 
scheme to general trajectories in the Schwarzschild space- 
time should be straightforward, based on the explicit ex- 
pressions given above for the Green's function, Eqs. (|]) 
and (^0|)-(|2^). In Ref. jl2| we further extend the regular- 
ization scheme to a more general class of static spherically 



** For circular orbits, F r (ALD) = 0. Thus, Eq. (||) is valid 
for the r component in circular motion as well. From the 
analysis of Ref. jlo| (concerning the evolution rate of E) it 
follows that in circular motion Eq. (32) holds also for the t 
component (with B t = 0; see also (14 ). 



symmetric spacetimes. We expect that the generalization 
to the electromagnetic self-force will be almost straight- 
forward. The important task of generalizing the method 
to the Kerr case and to the gravitational self-force seems 
more challenging; yet, we believe it should also be possi- 
ble. 

Finally, we note that a closely related regularization 
approach, which is also based on the multipole expansion, 
is currently being studied by Lousto p7| . 

We are grateful to Lior Burko for many interesting dis- 
cussions and helpful comments. One of us (A.O.) wishes 
to thank Kip Thorne for his warm hospitality during a 
visit in Caltcch, where this work was initiated. 
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